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Cosmological parameters deduced from the Planck measurements of anisotropies in the cosmic 
microwave background are at some tension with direct astronomical measurements of various pa¬ 
rameters at low redshifts. Very recently, it has been conjectured that this discrepancy can be 
reconciled if a certain fraction of dark matter is unstable and decays between recombination and 
the present epoch. Herein we show that if the superheavy relics have a branching into neutrinos 
~ 3x 10 -9 , then this scenario can also accommodate the recently discovered extraterrestrial 
flux of neutrinos, relaxing the tension between IceCube results and Fermi LAT data. The model is 
fully predictive and can be confronted with future IceCube data. We demonstrate that in 10 years 
of observation IceCube will be able to distinguish the mono-energetic signal from X decay at the 
3<J level. In a few years of data taking with the upgraded IceCube-Gen2 enough statistics will be 
gathered to elucidate the dark matter-neutrino connection at the 5cr level. 


We propose an explanation for the origin of Ice- 
Cube neutrinos |1] assuming heavy decaying dark matter, 
which is characterized by a lifetime and a comoving num¬ 
ber density that may resolve the tension between Planck 
data and low redshift astronomical measurements [ 2 |. 
The ensuing discussion will be framed in the context of 
dissipative dark matter. We consider a multicomponent 
hidden sector, with strong self interactions, featuring a 
massless dark photon which mixes with the ordinary pho¬ 
ton and a dark Z' which mixes with the Z. The dark 
photon induces a one-loop photon wave function mixing 
term eF^F'^ [3]. Early universe cosmology and galactic 
structure considerations constrain the dark photon mix¬ 
ing strength, e < 10 -9 [3]. Bounds on the Z — Z 1 mixing 
angle are less restrictive [5]. 

Cosmic Neutrinos — Recently, the IceCube Collab¬ 
oration reported the discovery of extraterrestrial neutri¬ 
nos p]. By establishing a strict veto protocol, the collab¬ 
oration was able to isolate 36 events in 3 years of data, 
with energies between 30 TeV < E v < 2 PeV. These 
events follow the expected spectral shape (oc E~ 2 ) of a 
Fermi engine, and are consistent with an isotropic distri¬ 
bution in the sky. A purely atmospheric explanation of 
the data can be excluded at 5.7 < 7 . 

At E]f s ~ 6.3 PeV, one expects to observe a dramatic 
increase in the event rate for v e in ice due to the Glashow 
resonance in which v e e~ —> W~ -A shower greatly in¬ 
creases the interaction cross section [ 6 ;. The hypothesis 
of an unbroken power law oc E~ a then requires a > 2.45 
to be consistent with data at ler J], More recently, the 


IceCube search technique was refined to extend the neu¬ 
trino sensitivity to lower energies E u > 10 TeV [5j. A fit 
to the resulting data, assuming a single unbroken power 
law and equal neutrino fluxes of all flavors, finds a softer 
spectrum 

-2.46±0.12 


x GeV 1 cm 2 sr 1 s 1 (1) 

and already mildly excludes the benchmark spectral in¬ 
dex a = 2. 

The neutrino flux in ([I]) is exceptionally high by astro¬ 
nomical standards, with a magnitude comparable to the 
Waxman-Bahcall bound [5]. A saturation of this bound 
can only be achieved within astrophysical environments 
where accelerator and target are essentially integrated. 
Potential candidate sources are discussed in [TO] . These 
powerful sources produce roughly equal numbers of 7r°, 
7 r + and tt~ in the proton-proton beam dump. The 7 r° ac¬ 
companying the parents of IceCube neutrinos decay 
into 7 -rays, which are only observed indirectly after prop¬ 
agation in the extragalactic radiation fields permeating 
the universe. These 7 -rays initiate inverse Compton cas¬ 
cades that degrade their energy below 1 TeV. The relative 
magnitudes of the diffuse 7 -ray flux detected by Fermi 
LAT m can then be used to constrain the spectral in¬ 
dex, assuming the 7 -rays produced by the 7r°’s accompa¬ 
nying the n ± escape the source. Figure [j] shows that only 
a relatively hard injection spectrum is consistent with the 
data. Indeed, if IceCube neutrinos are produced through 
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FIG. 1. The open symbols represent the total extragalactic 
7 -ray background for different foreground (FG) models as re¬ 
ported by the Fermi LAT Collaboration m- For details on 
the modeling of the diffuse Galactic foreground emission in 
the benchmark FG models A, B and C, see m- The cumu¬ 
lative intensity from resolved Fermi LAT sources at latitudes 
| 6 | > 20° is indicated by a (grey) band. The solid symbols 
indicate the neutrino flux (per flavor) reported by the Ice- 
Cube Collaboration. The blue points are from the data sam¬ 
ple of the most recent IceCube analysis [ 8 L The light grey 
data points are from the 3-year data sample of pQ, shifted 
slightly to the right for better visibility. The best fit to the 
data (extrapolated down to lower energies), is also shown 
for comparison [ 8 ]. The dashed line indicates the mono- 
energetic signal from dark matter decay. Note that a plot¬ 
ting of E 2 ® = EdF/(dQ dAdt din E) versus In E conserves 
the area under a spectrum even after processing the electro¬ 
magnetic cascade. Thus, the area of the 7 r° contribution to 
the diffuse 7 -ray spectrum (total diffuse 7 -ray flux provides 
an upper bound) implies the low energy cutoff (upper bound) 
to the origin of the neutrinos. 


pp collisions in optically thin extragalactic sources, the 
7 -rays expected to accompany the neutrinos saturate the 
Fermi LAT flux for a ss 2.2 |12|. The overall isotropy of 
the observed arrival directions and the fact that a PeV 
event arrives from outside the Galactic plane disfavor a 
Galactic origin. Moreover, for the Galactic hypothesis 
one must consider another important caveat, namely that 
the expected photon flux in the PeV range m has been 
elusive m- 

The difficulties so far encountered in modeling the pro¬ 
duction of IceCube neutrinos in astrophysical sources fu¬ 
eled the interest in particle physics inspired models. By 
far the most popular model in this category is the decay 
of a heavy massive (~ few PeV) relic that constitute (part 
of) the cold dark matter (CDM) in the universe [15] • The 
lack of events in the vicinity of the Glashow resonance, 
implies the spectrum should decrease significantly at the 
energy of a few PeV. Spectra from dark matter decays 
always exhibit a sharp cutoff determined by the particle 
mass. Furthermore, the 3 highest energy events appear 
to have identical energies, up to experimental uncertain¬ 
ties. A line in the neutrino spectrum would be a smoking 
gun signature for dark matter. If the heavy relic also de¬ 


cays into quarks and charged leptons, the mono-energetic 
neutrino line may be accompanied by a continuous spec¬ 
trum of lower-energy neutrinos, which can explain both 
the PeV events and some of the sub-PeV events. All 
of these considerations appear to be in agreement with 
the data [16] . Even much heavier relic particles, with 
masses well above a PeV, can generate the required neu¬ 
trino spectrum from their decays if their lifetime is much 
shorter than the present age of the universe [IT. The 
spectrum of neutrinos is modified by a combination of 
redshift and interactions with the background neutrinos, 
and the observed spectrum can have a cutoff just above 
1 PeV for a broad range of the relic particle masses. In 
this Letter we will reexamine the idea of a dark matter 
origin for IceCube events. 

Ho Measurements — Another, seemingly differ¬ 
ent, but perhaps closely related subject is the emerg¬ 
ing tension between direct astronomical measurements 
at low redshift and cosmological parameters deduced 
from temperature fluctuations in the cosmic mi¬ 
crowave background (CMB). Strictly speaking, the TT, 
TE, EE spectra recorded by the Planck spacecraft 
when combined with polarization maps (lowP) de¬ 
scribe the standard spatially-flat 6 -parameter ACDM 
model {flfc/i 2 , Ocdm^ 2 j O s , t, n s , A s } with high preci¬ 
sion: (i) baryon density, fl b h 2 = 0.02225 ± 0.00016; 
(ii) CDM density, Dcdm ^ 2 = 0.1198 ± 0.0015; (in) an¬ 
gular size of the sound horizon at recombination, 0 S = 
(1.04077 ± 0.00032) x 10 -2 ; (iv) Thomson scattering 
optical depth due to reionization, r = 0.079 ± 0.017; 
(v) scalar spectral index, n s = 0.9645±0.0049; (vi) power 
spectrum amplitude of adiabatic scalar perturbations, 
ln(10 10 A s ) = 3.094 ± 0.034 [18]. Planck data also con¬ 
strain the Hubble constant h = 0.6727± 0.0066, the dark 
energy density Da = 0.6844 ± 0.0091, the amplitude of 
initial density perturbations = 0.831 ± 0.013, and the 
mass density parameter D m = 0.3156 ± 0.0091 Q Un¬ 
expectedly, the H 0 inference from Planck observations 
deviates by more than 2.5cr from the previous interpre¬ 
tation of the Hubble Space Telescope (HST) data (based 
on over 600 cepheids in host galaxies and 8 samples of 
SNe la) which leads to h = 0.738 ± 0.024, including 
both statistical and systematic uncertainties m- A sep¬ 
arate study by the Carnegie Hubble program using mid- 
infrared calibration of the cepheid distance scale based 
on data from NASA’s Spitzer Space Telescope yields 
h = 0.743±0.021 (2Dj. Besides, the interpretation of grav¬ 
itational lensing time delay measurements of the system 
RXJ1131-1231 points to h = 0.787±g;gJ| (3J. 

The tension between the CMB based determination of 
the Hubble constant and the h value inferred from direct 
low redshift measurements is intriguing and deserves fur¬ 
ther attention. On the one hand, the underlying source 
of discrepancy could be some systematic uncertainty in 


1 Throughout we adopt the usual convention of writing the Hubble 
constant at the present day as Hq = 100 h km s —1 Mpc 1 . 
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the calibration [22]. On the other hand, it could trace 
a deficiency of the concordance model of cosmology. In 
the spirit of [23], it has been recently conjectured that 
Planck- inspired ACDM paradigm can be reconciled with 
HST measurements if a subdominant fraction fx of CDM 
is unstable and decays rather quickly with respect to 
the present Hubble time [2]. The width of the unsta¬ 
ble component r A is an independent parameter of the 
model. By forcing the X particles to decay after recom¬ 
bination Tx is bounded from above. Moreover, the X 
is assumed to decay (dominantly) into invisible mass¬ 
less particles of the hidden sector and hence does not 
produce too many photons. A joint fit to Planck , su¬ 
pernova, and HST data reveals that the base ACDM 
model, with r A = 0 , is outside the 2 er likelihood con¬ 
tours in the (T.y , fx) plane [2, 23]. The data instead 
favor 0.03 < fx < 0.12. The mean value and la 
error derived from a maximum likelihood analysis are 
h = 0.716 ± 0.020 [2] - Interestingly, within the same pa¬ 
rameter range the model could also alleviate the emerging 
tension with the cluster data. (See, however, 25].) For 
example, for f x — 0.10 and F A - — 2000 kms _ 1 Mpc _1 
the corresponding values of (l m ~ 0.25 and as — 0.80 [2] 
are marginally consistent with the 2 cr allowed contours 
by the Planck cluster mass scale :26| and the extended 
ROSAT-ESO Flux Limited X-ray Galaxy Cluster Survey 
(REFLEX II) [27]. For smaller values of fx and or T A 
the values of f l m and as move closer to the base ACDM 
model p]. Next, in line with our stated plan, we take 
fx — 0.05 and T A ~ 5000 kms _ 1 Mpc _1 (favored by 
lensing constraints [23]) as benchmarks and investigate 
what would be the CDM fraction required to decay into 
the visible sector to accommodate IceCube observations. 

Bump-Hunting — The two main parameters charac¬ 
terizing the X particle are its lifetime t a ~ 6 x 10 15 s and 
its mass mx, which is a free parameter. We assume the 
neutrino produced via X decay is mono-energetic, with 
energy e v = mx/ 2. The neutrino energy distribution 
from X decay is given by dN v /dE v = N v S(E V — e^), 
where N v is the neutrino multiplicity. We further as¬ 
sume the dominant decay mode into the visible sector, 
contributing to neutrino production, is X —> vv and so 
N v = 2. 

The evolution of the number density of neutrinos 
n v (E v ,z) produced at cosmological distances in the de¬ 
cay of X particles is governed by the Boltzmann conti¬ 
nuity equation, 

d[n„/(\ + zY] _ d[HE v n v /{l + zf] , ^ 

dt ~ dE v [Z) 

together with the Friedman-Lemaitre equations describ¬ 
ing the cosmic expansion rate H[z) as a function of the 
redshift 2 . This is given by H 2 (z) = Hq [fl m (l+ 2 ) 3 -|-DA]. 
The time-dependence of the red-shift can be expressed 
via dz = — dt (1 +z)H. We have found that for the consid¬ 
erations in the present work neutrino interactions on the 
cosmic neutrino background can be safely neglected [28] ■ 


In (§, 

Qu(E v ,t) = ^ , (3) 

t x dE u 

is the source term, nx(t) = Yx s(t) e~ t ^ Tx is the number 
density of A', Bx^pp is the neutrino branching fraction, 
s(t) is the entropy density, and 


Y x = 3.6 x 1(T 9 


(lx h 2 
mx / GeV 


(4) 


is the comoving number density at the CMB epoch. By 
solving ([ 2 ]) we obtain the (all flavor) neutrino flux at 
present epoch to, 

$(£?„) =-£-»„(£?„, 0) (5) 


c N v Y x s{to)„ e-Wr* 

4tt t x E v H{U) 


1 +z(t m )=e v /E v 


with s(to) — 2.9 x 10 3 cm -3 . 

Maximization of ([5]) yields the energy relation for the 
peak in the spectrum, 


ppeak _ 1 TO */ 2 

" ~ 2 1 + z(t x ) ’ 


( 6 ) 


which sets the mass of the X. Since z(tx) — 18 to ac¬ 
commodate the PeV peak in IceCube’s neutrino spec¬ 
trum we take mx — 76 PeV. Now, from Q we obtain 
Yx ~ 2.8 x 10 -19 . Finally, we normalize the cosmic neu¬ 
trino flux per flavor using ([!]) . The intensity of the mono- 
energetic signal at the peak is taken as 60% of the flux 
reported by the IceCube Collaboration, yielding a neu¬ 
trino branching fraction Bx^vp ~ 3x 10~ 9 into all three 
flavors. The width, an output of the Boltzmann equation, 
is shown in Fig. [l] It is evident that the mono-energetic 
neutrino spectrum is in good agreement with the data. 
In particular, the flux suppression at the Glashow reso¬ 
nance, 3>(£i)) es )/<lJ(EP eak ) ~ 0 . 011 , is consistent with data 
at la. 

The model is fully predictive and can also be con¬ 
fronted with Fermi LAT data. It is reasonable to as¬ 
sume that Bx^e+e~ ~ Bx~>i/ e D e ~ ^X—>uu ~ Bx~>dd m 
About 1 /3 of the energy deposited into either uu or dd is 
channeled into 7 -rays via 7r° decay and about 1/6 of the 
energy is channeled into electrons and positrons. As pre¬ 
viously noted, the 7 -rays, electrons, and positrons trigger 
an electromagnetic cascade on the CMB, which develops 
via repeated e + e” pair production and inverse Compton 
scattering. As a result of this cascade the energy gets re¬ 
cycled yielding a pile up of 7 -rays at GeV < E 1 < TeV, 
just below the threshold for further pair production on 
the diffuse optical backgrounds. We have seen that under 
very reasonable assumptions the energy deposited into 
neutrinos is comparable to the energy deposited into the 
electromagnetic cascade. Therefore, the neutrino energy 
density at the present epoch, 


CO v 


J E v n v {E v , 0 ) dE v = 2.6 x 1 CT 11 eV cm" 3 , (7) 
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TABLE I. Event rates (yr 4 ) at IceCube for E™ m = 1 PeV. 


E” ax /PeV 

spectrum oc E„ 2 46 

X decay spectrum 


Ve 

O* 

v T 

Ve 

O* 

v T 

2 

0.20 

0.18 

0.20 

0.25 

0.23 

0.25 

3 

0.27 

0.24 

0.27 

0.46 

0.41 

0.46 

4 

0.31 

0.27 

0.30 

0.58 

0.51 

0.56 

5 

0.34 

0.29 

0.32 

0.66 

0.56 

0.62 


provides reliable estimate of the cascade energy density 
(wcas ~ w„), which is bounded by Fermi LAT data to not 
exceed ~ 5.8 x 1CT 7 eVcm " 3 [29]. We conclude 

that the 7 -ray flux associated with the neutrino line is 
found to be about 4 orders of magnitude smaller than 
the observed flux in the Fermi LAT region. 

We now turn to discuss the possibility of distinguishing 
the neutrino line from an unbroken power-law spectrum 
without the neutrino line, with future IceCube data. The 
value of the spectral index is determined by the “low en¬ 
ergy” events. Following the best IceCube fit we adopt a 
spectrum a E~ 2A6 . We assume that the IceCube events 
below 1 PeV have an astrophysical origin. Indeed, the 
steep spectrum ex E~ 2A6 may suggest we are witnessing 
the cutoff of TeV neutrino sources running out of power. 
Using the IceCube aperture for the high-energy start¬ 
ing event (HESE) analysis [lj we compute the event rate 
per year above 1 PeV for both the neutrino flux given 
in ([I]) and that of |5]). The results are given in Table [I] 
As expected, the predictions from X decay are in good 
agreement with existing data. Because of the smeared 
energy-dependence of muon tracks, in what follows we 
will only consider cascades and double bang topologies 
initiated by charged current interactions of electron and 
tau neutrinos, as well as all neutral current interactions 
processes. We identify the events coming from the power 
law spectrum A/"n with background and adopt the stan¬ 
dard bump-hunting method to establish the statistical 
significance of the mono-energetic signal. To remain con¬ 
servative we define the noise = y/Afs + As, where A/s is 
the number of signal events. In 10 years of operation the 
total detection significance, 

5det = VA/B+A/-S ’ (8) 

would allow distinguishing the neutrino line from a sta¬ 


tistical fluctuation of a power law spectrum oc E~ 2A6 
at the 3 a level. Note that the shape of the distribu¬ 
tion with energy conveys additional information allowing 
one to distinguish the line signal from fluctuations of a 
power-law background. The proposed IceCube-Gen2 ex¬ 
tension plans to increase the effective volume of IceCube 
by about a factor of 10 [30] ■ This facility will not only in¬ 
crease the HESE sensitivity but also improve the energy 
resolution for muon tracks. In a few years of operation 
IceCube-Gen2 will collect enough statistics to elucidate 
the dark matter-neutrino connection with S ^ e t > 5cr. 

We end with an observation: IceCube data can also 
be fitted by a neutrino line peaking at E„ ~ 20 TeV 
superimposed over a power law spectrum (ex E~ 1 2 ) of as¬ 
trophysical neutrinos m ■ By duplicating our discussion 
for mx ~ 1 PeV it is straightforward to see that the 
model can also accommodate this neutrino line. 

Conclusions — We have shown that the PeV flux of 
extraterrestrial neutrinos recently reported by the Ice- 
Cube Collaboration can originate through the decay of 
heavy dark matter particles with a mass ~ 76 PeV and 
a lifetime ~ 6 x 10 15 s. On a separate track, the ten¬ 
sion between Planck data and low redshift astronomical 
measurements can be resolved if about 5% of the CDM 
component at CMB epoch is unstable. Assuming that 
such a fraction of quasi-stable relics is responsible for the 
IceCube flux we determined the neutrino branching frac¬ 
tion, Bx^vv ~ 3 x 10~ 9 . The model has no free param¬ 
eters and will be tested by future IceCube data. Indeed 
10 years of data taking will be required to distinguish 
the neutrino line from an unbroken power-law spectrum 
at the 3u level. The upgraded IceCube-Gen2 will collect 
enough statistics to elucidate the dark matter- neutrino 
connection at the 5cr discovery level in a few years of 
operation. 
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